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Quantum spin networks can be used to transport information between separated registers in a 
quantum information processor. To find a practical implementation, the strict requirements of ideal 
models for perfect state transfer need to be relaxed, allowing for complex coupling topologies and 
general initial states. Here we analyze transport in complex quantum spin networks in the maximally 
mixed state and derive explicit conditions that should be satisfied by propagators for perfect state 
transport. Using a description of the transport process as a quantum walk over the network, we 
show that it is necessary to phase correlate the transport processes occurring along all the possible 
paths in the network. We provide a Hamiltonian that achieves this correlation, and use it in a 
constructive method to derive engineered couplings for perfect transport in complicated network 
topologies. 
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I. INTRODUCTION 

In the quest toward a scalable quantum computer [1] , a 
promising model comprises distributed computing units 
connected by passive wires that transmit quantum in- 
formation [2-6]. This architecture would provide several 
advantages, since the wires require no or limited con- 
trol, easing the fabrication requirements and improving 
their isolation from the environment. For a simpler in- 
tegration in a solid-state architecture, the wires can be 
composed of spins. Following seminal work by Bose [7], 
which showed that spin chains enable transporting quan- 
tum states between the ends of the chain, the dynamics 
of quantum state transfer has been widely studied (see 
Ref. [8] for a review) and protocols for improving the 
fidelity by coupling engineering [9-16], dual-rail topolo- 
gies [17], active control on the chain spins [18] or on 
the end spins only [19-22] have been proposed. Re- 
cently these studies have been extended to mixed state 
spin chains [23-27], which are more easily obtained in 
high-temperature laboratory settings - making them im- 
portant protagonists in practical quantum computing. 
A further challenge to experimental implementation of 
quantum transport is the lack of chains with the desired 
coupling strengths, since coupling engineering is limited 
by fabrication constraints and by the presence of long- 
range interactions. These challenges highlight the need 
for a systematic study of mixed state transport in quan- 
tum systems beyond chains, including more complex net- 
work topologies. These topologies reflect more closely 
actual experimental conditions as well as systems occur- 
ring in nature. For example, there is remarkable recent 
evidence [28, 29] that coherent quantum transport may 
be the underlying reason for the high efficiency (of over 
99%) of photosynthetic energy transfer [30]. 
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To derive explicit conditions for perfect transport we 
quantify geometric constraints on the unitary propagator 
that drives transport in an arbitrary network. As a con- 
sequence, we find that transferring some mixed states in 
a network in general requires fewer conditions than pure 
state transfer. Transport conditions for pure states have 
been previously quantified [31], however, our method - 
relying on decomposing the propagator in orthogonal 
spaces - is fundamentally different and more suitable for 
mixed states. 

Perfect transport occurs when the bulk of the network 
acts like a lens to focus transport to its ends. To make 
this physical picture more concrete, we describe mixed 
state transport as a continuous quantum walk over the 
network [32-34] that progressively populates its nodes. 
Through this formalism, we derive constructive condi- 
tions on the Hamiltonian that results in the correlation 
of transport processes through different possible paths 
in the network. The correlation of transport processes 
leads to their constructive interference at the position of 
the two end-spins, giving perfect transport. While sim- 
ilar walk models have been applied to coherent transfer 
before (see [35] for a review), our work provides their first 
extension to transport involving mixed states. 

The insight gained by describing quantum transport as 
correlated quantum walks can be used to construct larger 
networks where perfect transport is possible. Here we 
show a strategy to achieve this goal by engineering the 
coupling strengths between different nodes of the net- 
work to construct weighted spin networks that support 
perfect mixed state transport. Feder [36] had considered 
a similar problem for pure states by mapping the quan- 
tum walk of N spinor bosons to a single particle; this 
has been extended in more recent work [37-39]. Here we 
find far more relaxed weighting requirements for mixed 
state transport, thanks in part to a fermionic instead of 
bosonic mapping. In turn, this could ease the fabrication 
requirements for coupling engineering. 

The paper is organized as follows. In Sec. II we define 
the problem of mixed state transport. Sec. Ill provides 
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the geometric conditions on the propagator for perfect 
transport in arbitrary networks. We finally present in 
Sec. IV the quantum walk formalism, which allows the 
correlation of transport processes over different paths and 
the construction of families of weighted networks that 
support perfect transfer. 



II. TRANSPORT IN MIXED-STATE 
NETWORKS 

Consider an iV-spin network A/", whose vertices (nodes) 
V represent spins and whose edges £ = {a^} describe 
the couplings between spins i and j (see Fig. f). The 
system dynamics is governed by the Hamiltonian H = 
J2i<j a ijHij, where H^ is the operator form of the in- 
teraction. In the most general case, a spin of Af may be 
coupled to several others, for instance, in a dipolar cou- 
pled network, ay ~ i s a function of the distance 
between the spins in the network. 

We assume that we can identify two nodes, labeled 1 
and N, that we can (partially) control and read out, in- 
dependently from the "bulk" of the network, and thus 
act as the " end" spins between which transport will oc- 
cur. The rest of the spins in the network can at most 
be manipulated by collective control. This also imposes 
restrictions on the network initialization [24-26]. To re- 
lax the requirements for the network preparation, we as- 
sume to work in the infinite-temperature limit [25] - a 
physical setting easily achievable for many experimen- 
tal systems - where the bulk spins are in the maxi- 
mally mixed state, pal. We will then consider the 
transport of a slight excess polarization from node 1 to 
node N. The initial state is p, ~ (1 + SZ%), where 
Z\ is the Pauli matrix acting on spin 1 and S <C 1 
denotes the polarization excess. Since only the trace- 
less part of the density matrix evolves in time, we will 
monitor the transport from pf" — Z\ to a desired final 
state p A = Zn- The fidelity of the transport process 

is then defined as F[t) = Tr{p A (t)Z N )/Tr{z\z 1 ), with 
p A (t) = U{t)pfW(t) being the evolved state. 

The polarization behaves like a wave-packet traveling 
over the network N [40, 41]. In most cases, the Hamilto- 
nian H drives a rapidly dispersive evolution, where the 
wave-packet quickly spreads out into many-body corre- 
lations among the nodes of A/", from which it cannot be 
recovered [34]. This is for example the case of evolu- 
tion under the naturally occurring dipolar Hamiltonian, 
which induces a fast-decay of the spin polarization as 
measured in solid-state NMR, even if many-body corre- 
lations can be detected at longer times [42]. 

In order to drive a dispersionless transport, thus en- 
suring perfect fidelity, the network Hamiltonian should 
satisfy very specific conditions. In this paper we will 
investigate these conditions by answering the questions: 
(i) What are the possible operator forms of the Hamilto- 
nian Hij for dispersionless transport? (ii) What are the 
coupling topologies and (iii) strengths that support 




FIG. 1: (Color online) Transport in a spin network. The 
network edges represent the interaction among spins (nodes) . 
The pink (dark) nodes represent the end spins between which 
transport should occur. The shaded region comprising yellow 
(light) nodes is the bulk network, Hb- An initial polarization 
packet (hatched) is prepared on spin 1 and allowed to propa- 
gate on the network through various possible paths (arrows) . 
For perfect transport, the polarization refocuses at spin N . 



perfect transport? 



III. CONDITIONS FOR PERFECT 
TRANSPORT 

A. Fidelity of mixed-state transport 

The condition for perfect transport, F = 1, can be ex- 
pressed in a compact form by using the product-operator 
(PO) basis [43]. For the iV-spin network system there are 
2 2N basis elements, 

B = Bi <K> Bbuik ® B w = {t,X u Y u Z 1 ,X 2 ,...,X 1 X 2 , 
. . . , ZiZ 2} XiX 2 X 3 , Z 1 Z 2 ■ ■ ■ Z N _iZ N }, (1) 

where B^n and Bbuik are the basis for the end and the 
bulk spins, respectively. 

Using the PO basis, the propagator U(t) = e~ lHt can 
be represented by a vector \U) = \cb x , cb 2 , ■ ■ ■ , cb 4N ] T 
in the 2 2N dimensional Hilbcrt-Schmidt (HS) operator 
space spanned by B [44]: 



U(t)=J2 i CB i {t) B h with c Bl 



Tr(BlBi) 



From an initial state with a polarization excess on spin 1, 



P~ 



Zi, the system evolves to 



P f = u P ?tf = Eij '•//'/.- />•,/,/>', , (2) 

yielding the transport fidelity to spin TV 

— n — T Tr Si i c* B .CB i B i Z 1 BjZ N =J2i c B,c* B ., 
with B 3 = ±Z 1 Z N B l , for [Bi, Z N ] T = 0. (3) 
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The last equation follows from the property that all ele- 
ments of B, except 1, are traceless. 



The fidelity derived in Eq. (3) has a simple form in the 
HS space. Note that in this operator space, the product 
Bill is a linear transformation, T : \U) — > Pb ( \U), where 
is a permutation matrix corresponding to the action 
of B>i [44] (here and in the following we denote operators 
in the HS space by a hat). The unitarity of U yields the 
conditions: 

\\U\\ = (U\U) = 1 (U\P Bt \U)=0, Bi^l. (4) 

Let us partition the HS space in two subspaces Q, Q, 
spanned by the basis G and G, 

G = Bi®Bbuik®{l, Z N }, G = BiOBbuiuOjXjv, Y N }, 

and note that all elements of G commute with Zn, while 
all elements of G anti-commute with Z^. We will label 
by superscripts Q and Q the projections of operators in 
these subspaces. Using this partition, we can simplify 
the expression for the fidelity of Eq. (3) to obtain 




PrW) 



FIG. 2: (Color online) Geometric interpretation of the condi- 
tion for maximum Z\ — > Zn transport fidelity. The unitary 
U is represented as a vector \U) (red thick arrow) in the HS 

space, with components \U S ) and in the subspaces Q 

(represented by an axis) and Q (represented by the shaded 
grey plane) . Pr \ U) (blue dashed arrow) is the reflection of 
\U) about the Q axis. Maximum fidelity occurs only when the 

Pz 1 z N causes a ■k- rotation of \U S 



F=(U\P Zi z n Pr\U) , 



(5) 



where Pr is a reflection about Q and Pz x z N is block- 
diagonal in the {G, G} basis (since Z\Z^ e Q): 



Pr, 



t G 
-I s 



P: 



Z\Zj\ 



p 



o 



o 

pQ 

Z\Zn 



(6) 



Rewriting the fidelity as the inner product between 
two vectors, F — {{Pz 1 z n U)\PrU) , provides a simple 
geometric interpretation of the perfect transport condi- 
tion, as shown in Fig. 2. The vector Pz ± z N \U) should 
be parallel to Pr | U) , which can be obtained if Pz 1 z N ro- 
tates by an angle n, while leaving \ U G ) unaffected. 

Alternatively, since Pz 1 z N just describes a tt— rotation of 
the vector \U) about the Z\Zm axis, for perfect transport 
the rotation-reflection operation S = Pz x z n Pr should be 
a symmetry operation for | U) . 

From Eq. (4) we have (U\Pz 1 z N \U) — and using 
Eq. (5) we can derive explicit conditions to be satisfied 
by the propagator to achieve perfect transport, F = 1 
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that simplify to 
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(7) 
(8) 

(9) 



When is this equation satisfied? By symmetry, it hap- 
pens when ||£/ e || = ||U^|| = 1/2, and \U G ) and 



(up to a phase) eigenvectors of P§ 1 z N an< ^ ^z x z 
eigenvalues ±1 respectively: 



U g ) are 
with 
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Z\Zjs 



\u G ) = + \u G ) 
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u L - 



(10) 



To enable perfect transport, \U) must thus have an equal 
projection on the two subspaces Q and G, as shown geo- 
metrically in Fig. 2. Also, intuitively from the symmetry 
operation S, all components of \U) lying on the plane Q 
should be rotationally symmetric with respect to Z\Zj^, 
while components of | U) lying on the plane Q should have 
reflection symmetry about Z\Z^. 

Note that Eq. (10) imposes fairly weak constraints on 
the transport unitaries, as opposed to the constraints for 
pure state transport [31, 45]. In particular Eq. (10) pro- 
vides no explicit constraint on the bulk of the network. 
For example, the two propagators, 

U\ = -Bbuik(l ± ZiZn) + B hulk (XiXjy ± YiYjsr), 

U 2 = B hulk (l±Z 1 Z N ) + B' hulk {X 1 Y N TY 1 X N ), [ ' 

with i?buik and B hvik arbitrary equi-norm operators (e 
span{Bbuik}) acting on the bulk, support perfect trans- 
port. Other propagators can be obtained thanks to an 
invariance property that we present in the next section. 
More generally, in Appendix A we explicitly provide a 
prescription to construct classes of unitaries for perfect 
mixed state transport. 



B. Invariance of transport Hamiltonians 

The fidelity F in Eq. (3) is invariant under a trans- 
formation U — VU, where V is unitary and commutes 
with S, that is, 



[V,P Zi z n Pr]=0. 



(12) 



This invariance can be used to construct Hamiltonians 
that support perfect transport starting from known ones. 
Consider an Hamiltonian H that generates the transport 
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evolution U — exp(-iHt). Then the transport driven 
by H is identical to that generated by the Hamiltonian 
H' = V^HV, where V satisfies Eq. (12). 

Ref. [46] proved similar symmetry requirements for 
Hamiltonians that transport pure states; here, however, 
we derived these Hamiltonian properties just from the 
geometric conditions on U. Ref. [47] treated a similar 
problem, defining classes of Hamiltonians that perform 
the same action on a state of interest. A special case of 
this result was used in [23] to study transport in a mixed- 
state spin chain driven either by the nearest-neighbor 
coupling isotropic XY Hamiltonian, 

Hxy = ^T+ +1 with T+ = (S+S- + Sr S +) (13) 

or double-quantum (DQ) Hamiltonian 

#dq = V with D+ = (S+S+ + Sr S r), (14) 

with — \{Xj ± iYj). The unitary operator relating 
the two Hamiltonians, V = Y\ k > X k > , where the product 
k extends over all even or odd spins, does indeed satisfy 
Eq. (12). 



Quantum information transport via mixed state 
networks 



The requirements for perfect transport (Eq. (10)) can 
be easily generalized to the transport between any two 
elements of B, say from X to J 7 . One has simply to 
appropriately construct the subspaces Q and Q and the 
corresponding permutation operator Pxt- 

One could further consider under which conditions this 
transport (for example, X\ — > Xn) can occur simultane- 
ously with the Z\ — > Zj\f transport already considered. 
More generally, the simultaneous transfer of operators 
forming a basis for Bx would enable the transport of 
quantum information [10, 25] via a mixed-state network. 
The unitary U should now not only be symmetric un- 
der S, but should also under a similar operator derived 
for XiXpf. The requirements on U thus become more 
stringent and only a special case of the propagators con- 
structed in Eq. (A3) in Appendix A is allowed, 



U = Sbuik(l ± Z 1 Z N + X,X N ± Y 1 Y J 



N) 



(15) 



This is exactly a SWAP operation (up to a phase) be- 
tween the end-spins, which can also lead to a transfer 
of arbitrary pure states between 1 and N. Therefore, 
we find that perfect transport of non-commuting mixed 
states between the end-spins also allows transport in 
pure state networks. We note that quantum information 
could be encoded in multi-spin states [25, 48] that sat- 
isfy proper symmetry conditions and thus do not impose 
additional conditions on the transport propagators. 



D. Which Hamiltonians support mixed state 
transport? 

It would be interesting to determine which Hamiltoni- 
ans can generate propagators \U(t)) = exp(-iHt) |U) for 
perfect transport. Unfortunately, deriving requirements 
for the Hamiltonian from the conditions on the unitaries 
is non-trivial; however, as we show below, one can still 
extract useful information. 

A general Hamiltonian can be decomposed as H = 
H s + , where H S,G lie in the subspaces Q and Q, re- 
spectively. We cannot set H = since the Hamiltonian 
does need to have a component that is non-commuting 
with the target operator (Zn) in order to drive the trans- 
port. If H = H s , odd powers of H are in Q, while even 
powers of H belong to Q. Then the propagator has con- 
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H \H) 
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We can demonstrate that in this case the Hamiltonian 
must satisfy two conditions to drive perfect transport. 
First, the "vector" form of the Hamiltonian must be an 
eigenstate of Pz x z N , Pz t z N \H) = — \H), which ensures 
that the second equation in (10) is trivially satisfied, as 

Pi 



Z-i Z A 



U^). Second, since we have 



P Zl z N \U G ) - \Z X Z 



N 



2! 



H\H) 



(jty 

4! 



H 6 \H) 



the first equation in (10) implies that H 2n = \ (\—Z\Zn) 
for any n. 

These conditions are for example satisfied by the XY- 
likc Hamiltonian, H = -BbuikT-jjy, where -Bbuik is any op- 
erator acting on the bulk and T^ N = (S^ ± S^S^). 
In this case, at t — ir/4, all conditions in Eq. (10) are sat- 
isfied and perfect transport is achieved. Indeed the XY 
Hamiltonian has been widely studied for quantum trans- 
port [7, 9] and it is interesting that we could derive its 
transport properties solely by the symmetry conditions 
on the propagator. 

An Hamiltonian H = HY with support only in Q is 
however a very restrictive case as it refers to the situa- 
tion where all nodes of the network are connected to N. 
Hamiltonians with support in both subspaces are more 
experimentally relevant, as they correspond to a common 
physical situation, where the ends of the network are sep- 
arated in space and direct interaction between them is 
zero or too weak. In the following, we will consider this 
more general situation, although restricting the study to 
XY Hamiltonians in order to derive conditions for perfect 
transport. 
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IV. PERFECT TRANSPORT IN NETWORKS: 
CORRELATED QUANTUM WALKS 

In the following, we will consider the network Af to 
consist of spins that are coupled by XY-like interac- 
tion, {T^j}. We focus on this interaction since it has 
been shown that with appropriate engineered coupling 
strengths, <Xij oc Wi(N — i)5j i+i, the XY-Hamiltonian 
can support perfect transport in linear spin chains (see 
e.g. [10-12, 49]). Thanks to the invariance property 
described in Sec. IIIB, this analysis applies to a much 
broader class of Hamiltonians, in particular to the DQ 
Hamiltonian. 

We assume that the end spins of J\f are not directly 
coupled, thus transport needs to be mediated by the bulk 
of the network. The simplest such topology is a A-type 
configuration where the end spins are coupled to a single 
spin in the bulk. The Hamiltonian Aj = {T^+T+ N ) /y/2, 
where j is a spin in the bulk, is enough to drive this 
transport. In this case, (A|) n = A| Vn, and hence the 
propagator is 



U = exp(-iHt) = 1 + [cos(t) 



l]H 2 



ism(t)H, (17) 



where H 2 = 1/2[(T+) 2 + {T+ N f + T+ N ). Since (1-2H 2 ) 
has the form of U3 in Eq. (A3), setting t — it ensures 
U = U3, yielding perfect transport. This is an expected 
result, since this simple lambda- network is just a 3-spin 
linear chain. This result can be extended to longer chains, 
as long as engineered couplings ensure that the resulting 
Hamiltonian is mirror-symmetric [10, 16, 45]. 

A different situation arises when there is more than one 
transport path possible, that is, the end spins are coupled 
to more than one spin in the bulk with an Hamiltonian 
H = X^'ebuik a j^j- For example, Fig. 3 depicts a net- 
work similar to the one considered in [26, 50] where there 
are three A paths between the end-spins. Even if each 
path individually supports perfect transport, evolution 
along different paths may not be correlated, leading to 
destructive interference reducing the fidelity (see Fig. 3). 

Perfect fidelity can be achieved only if different paths 
can be collapsed into a single "effective" one that sup- 
ports perfect transport (Fig. 8). This strategy not only 
allows us to determine if an Hamiltonian can support per- 
fect transport, but it also gives a recipe to build allowed 
Hamiltonians, by combining simpler networks known to 
support perfect transport into more complex ones. 

To this end, we use the fact that linear chains enable 
perfect transport with appropriate engineered couplings. 
Our first step will then to give conditions under which 
two chains of the same length (with end-spins in com- 
mon) can be combined. To obtain these conditions, we 
describe the evolution of the spin polarization as a quan- 
tum walk over the operators in the network [33, 34]. This 
description reveals the need to correlate the parallel paths 
over the network, in order to achieve a constructive re- 
focusing of the polarization at the other end of the net- 
work. We then generalize the conditions by a recursive 
construction to quite general networks. 
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FIG. 3: (Color online) (a) A-type network with three A paths 
between the bulk and end-spins and equal coupling strength 
a. (b) Transport fidelity as a function of normalized time 
for the A-network coupled by the XY-Hamiltonian X£T (blue 

dashed) and the modified XY-Hamiltonian TjJ (red solid). 
In the latter case, correlated quantum walks lead to perfect 
transport. 



A. Transport as a quantum walk over J\f 

We describe the transport evolution as a quantum walk 
over the network, which progressively populates opera- 
tors in the HS space. This process of progressively pop- 
ulating different parts of the HS space upon continuous 
time evolution under the Hamiltonian can be considered 
as a quantum walk [33, 34, 51]. We first expand the 
transport fidelity F(t) (Eq. 5) in a time series, 



F(t) 



U 
1H 2 



Pz x z n Pr 



2! 



U 



H, 



it(u [h,P ZiZn P r ] 



U 



u 

(18) 



H, PzxZ n Pr 

with \Uq) = Defining the nested commutators, 

Co = P Zi z n Pr ; C n = [H,C n -x] , (19) 
Eq. (18) takes the form 



(20) 



n=0 



where the expectation value is taken with respect to \Uq). 

A large part of the Hamiltonian commutes with 
Pz ± z n Pr and can be neglected. We can isolate the non- 
commuting part by defining the operator A via the rela- 
tionship 



[H,Pz 1 z n Pr]=AP ZiZn P b 



(21) 



The operator A and its nested commutators C r f = 
[H,C^_i\ (with Cq — A) have a simple graphical con- 
struction. The commutation relations, 



ij ' jk 



-ZjTT ; [T+,T+]=0, 



(22) 



() 




FIG. 4: (Color online) Graphical representation of commu- 
tators for (a) the XY and (b) the modified XY Hamiltonian. 
The commutator between two top legs (blue) of the directed 
graph is the third edge (red). In case of the XY Hamiltonian, 
the commutator is conditioned on node j. 

(see Fig. 4) can be used to provide a simple prescription 
to graphically determine the flip-flop terms in C£. For 
any two edges, one in C^_ 1 and one in H, that share a 
common node, contains the edge required to complete 
the triangle between them. Thus, each higher order in 
the commutation expansion creates a link between nodes 
in the network, progressively populating it. We will refer 
to the operators as quantum walk operators, since as 
we show below, the nested commutators C n in Eq. (20) 
can be built exclusively out of them. 

Consider the network of Fig. 5(a), with coupling 
strengths ctij = 1: A contains only the edges of N that 
connect to node 1, as represented by the red lines in 
Fig. 5(b). Fig. 5(c-d) represent the higher order Commu- 
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FIG. 5: (Color online) (a) A six spin network with two paths 
Pi and P2 between the end spins, (b)-(f) represent graphically 
the successive orders of the quantum walk operators C„: A 
red line linking two nodes indicates that there is a flip-flop 
term JJy between them, while path-dependent prefactors are 
not depicted. Once the walk has covered the entire network, 
successive orders in reproduce C3 1 and C£ . The explicit 
expressions for the commutators are shown in Table I. 



--EE- E V)( V) 



For a commutator C n to yield a non-zero contribution 
to the fidelity, the product of the operators should 
be proportional to Z\Zjs[, that is, it should evaluate to 



tators, with a red line linking two nodes denoting a term 
Tf- between them. We note that the graphical construc- 
tion only predicts the presence of a flip-flop term link- 
ing two nodes in the commutator C£, while the explicit 
forms of the commutators is generally more complex, as 
shown in Table I, with additional appropriate weights for 
arbitrary coupling strengths a^ . Still, as we now show, 
only the T^j terms are important to determine the fi- 
delity, and the presence of a T 1N term in the graphical 
series is an indication that transport can occur between 
the end-nodes. 



The commutators C n can indeed be written in terms 
of the C£ nested commutators, 



C « = E( k ftx-kCk, (23) 



yielding an expression for the fidelity containing only 
products of the nested commutators C^: 



(24) 
I 

even powers of T^ N . Hence very few terms appearing in 
Eq. (24) actually contribute to the transfer fidelity F. 

The geometric construction of only yields the XY 
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XY Hamiltonian 


Modified XY Hamiltonian 
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5T+ + 5T+ + lST+v 


c 4 A 


-(61 + ZZ 1 Z 4 Z 5 Z N )T+ - {iSZ 1 Z 4 + 3Z 3 Z N )T+ - {6Z 1 Z 2 + 3Z S Z N )T+ 


-9T+ - 9T+ - 9T 3 + 




+2(1 + Z 1 Z 2 Z i Z 5 )T+ - (61 + 3Z 1 Z 2 Z 3 Z N )T+ + 2(1 + Z 1 Z 2 Z 3 Z i )T+ 


+4T 3 + - 9T+ + 4T+ 



TABLE I: Nested commutators C„ (walk operators) corresponding to the graphs in Fig. 5 if the edges represent the XY 
Hamiltonian or the modified XY Hamiltonian. In the first case, note the presence of path dependent Zj prefactors, which are 
absent if the modified XY Hamiltonian is used. This allows for the correlation of transport through parallel paths. 



operators contained in each commutator, but it does not 
reflect the appearance of prefactors oc Zj (due to the 
commutator in Eq. (22)) that are explicitly written out 
in Table I. Thus, the geometric construction gives a nec- 
essary condition for transport, but not a sufficient one. 

The operators describe the sum of walks over dif- 
ferent paths: for example in C 2 , Z 2 Z^T^ N can be in- 
terpreted as the information packet reaching node N 
through path Pi in Fig. 5(a), while Z^Z$T^ N represents 
propagation via path P 2 . These two terms could in prin- 
ciple contribute to the fidelity, as they contain T^ N . How- 
ever, the additional path dependent factors \\ k Z^ lead 
to a loss of fidelity. Transport through different paths 
yield different \\ k Z k factors, resulting into a destructive 
"interference" effect. Note also that since different paths 
are weighted by different correlation factors, they cannot 
be canceled through some external control to recover the 
fidelity. In the following section we show how a modified 
Hamiltonian can remove this path-conditioning and thus 
drive perfect transport. Note that the path-dependent 
factors are as well unimportant in the case of pure states, 
provided the states reside in the same excitation mani- 
fold [25, 52, 53]. 



B. Correlating quantum walks over JV": Modified 
XY-Hamiltonian 



To remove the path-conditioning one should modify 
the Hamiltonian so that the Zj term in the commutator 
Eq. (22) disappears. 

This can be done via a modified XY-Hamiltonian 



n z « 

i<u<j 

since it satisfies this condition: 



rji-\- T^i 
1 ij^ 1 jk 



^ik ' 



± ij ' ± kl 







(25) 



(26) 



that distinguishes paths between the two nodes i and j. 
Note that when the network J\f is a simple linear chain 
with nearest-neighbor couplings the modified Hamilto- 
nian Ty is equivalent to the bare XY Hamiltonian. The 
modification in Eq. (25) of the XY-Hamiltonian could 
also be seen as mapping the spin system into a set of non- 
interacting fcrmions [54, 55] via a Jordan- Wigner trans- 
formation [56, 57], since Ci — Y\ u<i Z u Sf are operators 
that satisfy the usual fermionic anti-commutation rela- 
tionships. When these modified operators are employed 
in the network Af of Fig. 5(a), the two paths Pi and P 2 
in Fig. 5(a) are indistinguishable or, equivalently, they 
become perfectly correlated (see Table I). In effect, the 
modified XY Hamiltonian drives the quantum walks over 
different paths through a common set of operators of B. 
This is shown in Fig. 6 for a simple A-network consisting 



(a) 

Walk 1 -> 2 -> N 



(b) 



Zi iT 12 Z 2 jiT 2N Z /N z 1 if- Z 2 iT - N Zj, 
-iZ 2 T{ N * 




These operators now depend on the number of nodes be- 
tween % and j, thus introducing a metric in the spin-space 



- iZ 3 T{ 
Z\ iT{^Z 3 iT 3 ~ N Z N 
Walk 1 -> 3 -> N 



FIG. 6: (Color online) Operators appearing in the quantum 
walk of a network consisting of two A paths, 1 — > 2 — > N 
and 1 — ► 3 — >■ N. In panel (a), where the transport is driven 
by the XY-Hamiltonian, the two paths through spins 2 and 
3 are different, as they traverse a different set of operators, 
and are thus depicted in two separated grey panels. In panel 
(b), where we consider the modified XY Hamiltonian, both 
walks go through a common set of operators. The previously 
separate walks are bridged by the operators in the red box, 
making both walks indistinguishable and hence correlated. 




FIG. 7: (Color online) Engineered spin networks, (a) The network of Fig. 5(a), correlated by means of the modified flip-flop 
Hamiltonian, can be engineered to yield perfect transport by weighting the coupling strength a with the coefficient shown. 
Here 7 is any positive parameter, 7 < f. The two paths Pi and P2 can be collapsed to form a linear chain, since nodes (2,4) 
and (3,5) are equivalent, (b) A more complicated network consisting of three paths Pi — P3 (black arrows), with engineered 
strengths parametrized by < 71,72 < 1. The three paths can be collapsed into an effective 4-spin linear chain, with equivalent 
nodes (2,5) and (3,4,6). (c) Transport fidelity as a function of normalized time for the network in (b) with 71 = 0.8,72 = 0.5, 
in case the edges are the XY {T^} or modified XY {J^t } operators. In this last case perfect fidelity is achieved, while for the 
usual XY-Hamiltonian the path-conditioned interference leads to poor transport fidelity. 



of two A paths. 

The graphic construction used to calculate the trans- 
port over the network in Fig. 5 remains unchanged, ex- 
cept that now the red lines between two nodes denote 
modified flip-flops between them. Crucially there are 
no path dependent prefactors and symmetric nodes in 
each path become equivalent in each of the operators C„. 
It is then possible to collapse different paths into a single 
effective one, until a complex network J\f is collapsed into 
a linear chain. This is depicted in Fig. 7(a). 

We can express this result more formally, by defining 
collapsed XY operators, where we denote in parenthesis 
equivalent nodes in two parallel paths: 



T: 



1 



7ik 



0^ 



(27) 



where 7^ and 7^ are arbitrary parameters, < 7^ , 7^ < 
1 (see also Appendix B). Remarkably, these operators 
satisfy the same path-independent commutation rela- 
tions as in Eq. (26) 



T. 



T, 



± 



(28) 



thus showing that intermediate equivalent nodes can be 
neglected in higher order commutators. In addition the 
nested commutators C^, and the graphical method to 
construct them (Fig. 4), remain invariant when substi- 
tuting the modified XY operator with the collapsed op- 



erators T. 



Using the collapsed operators, the network of Fig. 5(a) 
can thus be reduced to a simpler linear chain (Fig. 7). 
Analogous arguments for path-collapsing were presented 
in [58], and have been applied before to some classes 
of graphs [9, 38]. In the following we show that path- 
equivalence could be constructed even for more complex 
network topologies, since, as we described, path collaps- 
ing can be derived just from the commutation relation- 
ships between the edges of the network. 



C. Engineered spin networks 



The path collapsing described in the previous section 
provides a constructive way to build networks, with ap- 
propriate coupling geometries and strengths, that achieve 
perfect transport. Alternatively, given a certain network 
geometry, the method determines all the possible cou- 
pling strength distributions that leave its transport fi- 
delity unchanged. 

For example, starting from a linear chain, any node can 
be substituted by two equivalent nodes, thus giving rise 
to two equivalent paths. Then, within the subspace of the 
equivalent nodes, the couplings can be set using Eq. (27) 
with arbitrary weights 7, thus giving much flexibility in 
the final allowed network. The engineered network cor- 
responding to Fig. 5(a) is represented in Fig. 7(a), where 
equivalent nodes from P\ are weighted by 7, while those 



from P 2 are weighted by \f\ 



A more complex example is shown in Fig. 7(b), where 
the network is built combining the networks in Fig. 3(a) 
and Fig. 5(a). It consists of three paths and can be col- 
lapsed into a 4-spin linear chain. The couplings shown 
lead to perfect Z\ — > Zjy transport for arbitrary path 
weights 71 and 72, with < 71,72 < 1, as shown in 
Fig. 7(c). The network engineering scheme can be re- 
cursively integrated to construct larger and more com- 
plicated network topologies (see for example Fig. 8). 

Similar weighted networks have been considered be- 
fore for bosons [36] . The engineered couplings derived by 
mapping quantum walks of N spinor bosons to the walk 
of a single particle are however much more restricted than 
what we found here via the mapping of spins to non- 
interacting fermions. 
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FIG. 8: (Color online) A complex network topology that can 
be engineered for perfect transport with the modified XY 
Hamiltonian. Collapsing the equivalent nodes along the lines 
^2 — (-5 leads to three equivalent paths that can be suitably 
engineered for perfect transport. 



V. CONCLUSIONS AND OUTLOOK 

Experimental implementation of quantum information 
transport requires relaxing many of the assumptions 
made in ideal schemes. In this paper we analyzed a 
physical situation that is closer to experimental settings 
- information transport in mixed-state spin networks 
with complex topologies. We first derived general con- 
ditions on propagators that allow perfect transport in 
these mixed-state spin networks. We used the conditions 
on the propagators to show that there exist classes of 
symmetry transformations on the Hamiltonians driving 
the transport for which the transport fidelity is invari- 
ant. We also showed that the propagator conditions 
also imply that transporting some mixed states requires 
fewer control requirements than pure state transport, an 
added advantage to using mixed-state channels in quan- 
tum information architectures. 

In order to study quantum transfer in complex spin 
networks, we described the dynamics as a continuous 
quantum walk over the possible paths offered by the 
network. This description provided a graphical con- 
struction to predict the system evolution, which high- 
lighted the need of correlating the transport processes 
occurring along different paths of the network to ob- 
tain perfect transport. We thus introduced a modified 
XY-Hamiltonian, based on Jordan- Wigner fermioniza- 
tion, that achieves correlation among paths by establish- 
ing a metric for the quantum walks occurring on the net- 
work. Conversely, the graphical construction could be 
as well used to study the generation from the usual XY- 
Hamiltonian of states of interest in measurement-based 
quantum computation [59]. 

Finally, the quantum-walk picture and the graphical 
construction lead us to define a constructive method to 
build complex networks from simpler ones, with appro- 
priate coupling geometries and strengths, that achieve 
perfect transport. We thus found that there is consid- 
erable freedom in the choice of topology and interaction 



strength that still allows perfect transport in complex 
networks. While the requirement of a well-defined net- 
work topology could be further relaxed [60], the precise 
construction proposed in this paper would provide faster 
transport and the freedom in the coupling distributions 
could make these networks implementable in experimen- 
tal systems. 

Acknowledgment 

This work was partially funded by NSF under grant 
DMG-1005926. 

Appendix A: Constructing perfect transport 
unitaries 

Here we show how the conditions specified in Eq. (10) 
could be used to construct perfect transport propagators. 
Our motivation for this is to demonstrate that the con- 
ditions of Eq. (10) are very weak, in the sense that is 
possible to construct an infinite classes of unitaries that 
support Zi — > Zjq transport. 

Consider the matrix forms of PS 7 and P% 7 in the 

two-dimensional {l,iV} subspace of Q and Q: 

G ~ {{1, Z X Z N }, {Z u Z n }, {X x ,Y x Z n }, {Y U X X Z N }} 

g ~ {{XxXjy, Y 1 Y N }, {X 1 Y If ,Y 1 X I f}, {X Nl Z X Y N }, 

{Y N ,Z 1 X tf }} (Al) 

where the ~ refers to the restriction in the {l,iV} sub- 
space. Then, for this order of basis, the matrix forms are 
block diagonal 

H iZn = diag([X, X,-Y,Y}) 

P§ iZn = diag([-X, X,-Y,Y]) (A2) 

where X and Y are the standard Pauli matrices, whose 
eigenvectors with eigenvalues ±1 are respectively [1, ±1] T 
and [l,±i] T ; this imposes a restriction on U. If 
-Bbuik, -E"buik spanlBbuik}, one can explicitly list from 
Eq. (10) possible forms of U for perfect transport, 

Ui = B bulk(l ± Z\Zm) + B hulk (XiX]y ± YiYn), 

U 2 = B hulk (t±Z 1 Z N ) + B' hulk (X 1 Y N TY 1 X N ), 

U 3 = B hulk (Z 1 ±tZ N ) + B' hulk (X 1 X N ±Y 1 Y N ), 

U± = Bbuik(Zi ± Zn) + B bulk (XiYjv T YiXn), 

U 5 = B hul ±{X 1 ±iY 1 Z N ) + B' hulk (X N TiZ 1 Y N ), 

U 6 = B hvXk {X 1 ±iY 1 Z N ) + B' hu]k (Y N ±iZ 1 X N ), 

U 7 = B hulk (Y 1 ±iX 1 Z N ) + B' hulk (X N ±iZ 1 Y N ), 

U 8 = B baSk (Y l ±iXiZ N ) + Bl ulk (Y N TiZiX N ), 

(A3) 
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Note that the bulk of the network specified by -Bbuik and 
£? bulk can be any arbitrary operators with equal norms. 
In fact, the invariance described in Sec. IIIB could be 
used to show that the eight forms of U in Eq. (A3) are 
equivalent to U\ or XJi. 

Of course, one can combine the forms in Eq. (A3) 
to form other propagators that continue to support 
perfect transport. Consider for example a propaga- 
tor constructed out of XJ\ and U% in Eq. (A3), with 

Sbulk, 5 bulk = 1 



We have then the following product relations: 



(T±\ 2 _ r+ T ± T + _ 1 ( T ± 
\ ij) ~ ij ' ij jk ~ 2 V 'ifc 

and the commutation relations: 



ij ij ij' 



ZjTg) 



T ?j,T+ k 



+ 

ij ' 



T+, Zt\ =-2Tr jt [T^,Z t ]=-2T, 



(B4) 



(B5) 



U = X 1 (l±Z 1 Z N )+X 2 (X 1 X N ±Y 1 Y N )+X 3 (X 1 Y NT Y 1 X N ) 

where Xj are coefficients to be determined. Then, from 
Eq. (4) we have 



(U\U) 



1 



|Ai| 2 + |A 2 | 2 + |A 3 | 2 



1 



(U\P ZiZn \U) = => |Ai] 2 = ]A 2 ] 2 + |A 3 | 2 
(U\ P Zn \U) = Im(A^A 3 ) = (A4) 

Other conditions in Eq. (4) are satisifed trivially. 
Eq. (A4) can be solved exactly; for example Xj = 
{1/V2, 1/2, 1/2} is a solution. Importantly however, if 
the B bulk 's were different from each other for U\ and U2, 
the set of equations Eq. (A4) becomes far simpler. 

In summary, achieving Z\ — » Zn transport requires 
weak conditions on the propagator driving the transport. 
This is as opposed to perfect pure state transport, that 
requires the propagators to be isomorphic to permutation 
operators [31] that are mirror symmetric [45] about the 
end spins of the network. 



Appendix B: Properties of flip-flop and 
double-quantum Hamiltonians 

In this appendix, we present simple relations satisfied 
by the flip-flop (XY) operators that will be used in the 
main paper. Note that the double-quantum (DQ) oper- 
ators in Eq. (14) follow analogous equations. In what 
follows, distinct indices label distinct positions on the 
spin network unless otherwise specified. We start with 
the definition of the operators S and E: 



1 



1 



E t = 2 (1 ± Z A ' S t = 2 {Xl ± iYj) - (B1) 
These operators satisfy the following product rules: 



Z j Sf=±Sf, (Sf? = EfEj = 0, 
SfSf = (Eff = Ef. 



(B2) 



We define the flip-flop operators X~ and Lfy 



T: 



{E+Ej±Er E +) (B3) 



From the definition in Eq. (B3) it follows that 



T±=±T±- Zj T+=T h 



We define the modified flip-flop operators T^, 



Tf 



n z *> 



(B6) 



i<u<j 



obtained by multiplying the flip-flop operator in Eq. (B3) 
by a factor of Z u for all nodes between i and j. The 
modified flip-flop operators follow especially simple com- 
mutation rules 



rp-\- T^i 

ij ' jk 



- -T T ■ 

~ ± ik ' 



rp-\- rp± 



= 



(B7) 



Note that crucially, these commutators only depend on 
the initial and final nodes (i and fc), and are indepen- 
dent of intermediate nodes. In a physical analogy, the 
modified operators behave as if they were path inde- 
pendent. Thus, when considering two (or more) paths, 
we could omit any intermediate node, since it would 
not enter in the ensuing commutators. We then denote 
equivalent nodes in parenthesis -for example, (j, k) means 
nodes j and A; are equivalent- and define the collapsed 
operators: 



T± = 

i{j,k) 



1 



7^+7^ 



\lij T ij + 7»fc^fc 



(B8) 



where 7y, 7^ are arbitrary parameters, 0<jij, r fik< 1. 
The collapsed operators satisfy commutation relations 
similar to Eq. (B7): 



± iU*)> U,k) 



1 a 1 



T, 



{j,k)V it 



T+ 



i{r, 



(j,k)(m,n) 



(B9) 



The collapsed operators in Eq. (B8) can be generalized. If 
I = (01, 02, • • ■ , a m ) and J = (61, 62, • • ■ , b n ) denote two 
sets of equivalent nodes, we have the collapsed operator 

1 



T ± _ 



:EEw^' ( B1 °) 



Z^i=l L~ij = l l ai bj 1=1 J = l 

which satisfies the commutation relationships: 

= 0. 



IJ ' JK 



rp-\- rp± 

IJ ' KL 



(Bll) 
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